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Abstract 

We present a new perturbative formulation of non-equilibrium thermal field theory, based upon non-homogeneous free propagators 
C^ and time-dependent vertices. The resulting time-dependent diagrammatic perturbation series are free of pinch singularities without 
'""' the need for quasi-particle approximation or effective resummation of finite widths. After arriving at a physically meaningful 



(N 



definition of particle number densities, we derive master time evolution equations for statistical distribution functions, which are 
valid to all orders in perturbation theory and to all orders in a gradient expansion. For a scalar model, we perform a loopwise 
Q^truncation of these evolution equations, whilst still capturing fast transient behaviour, which is found to be dominated by energy- 
.^ violating processes, leading to the non-Markovian evolution of memory effects. 

\^ Keywords: non-equilibrium thermal field theory, non-homogeneous backgrounds 
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1. Introduction 

The description of out-of-equilibrium many-body field- 
theoretic systems is of increasing relevance in theoretical and 
experimental physics at the density frontier. Examples range 
from the early Universe to the deconfined phase of QCD, the 
quark-gluon plasma, relevant at heavy-ion colliders, such as 
RHIC and the LHC as well as the internal dynamics of com- 
pact astro-physical phenomena and condensed matter systems. 

In this Letter, we present the key concepts of a new per- 
turbative approach to non-equilibrium thermal quantum field 
theory, where master time evolution equations for macroscopic 
observables are derived from first principles. A comprehen- 
sive exposition of this new formulation is provided in |[T|. In 
contrast to semi-classical approaches based on the Boltzmann 
equation |2-9|, this new approach allows the systematic incor- 
poration of finite-width and off-shell effects without the need 
for effective resummations. Furthermore, having a well-defined 
underlying perturbation theory that is free of pinch singulari- 
ties, these time evolution equations may be truncated in a per- 
turbative loopwise sense, whilst retaining all orders of the time 
behaviour Several studies appeared in the literature 1 10 - 401 
proposing quantum-corrected transport equations, based upon 
systems of Kadanoff-Baym equations |'4T1, functional renor- 
malization group approaches 1421 or expansion of the Liouville- 
von-Neumann equation P3l l44l . Whilst retaining all orders in 
perturbation theory, the existing approaches often rely on the 
truncation of gradient expansions 145, 461 in time derivatives, 
quasi-particle approximations or ad hoc ansaetze in order to ob- 
tain calculable expressions or extract meaningful observables. 
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In this new perturbative formalism, the loopwise-truncated evo- 
lution equations are built from non-homogeneous free propaga- 
tors and time-dependent vertices. This diagrammatic approach 
encodes both spatial and temporal inhomogeneity already from 
tree-level, without resorting to any such approximations. 

2. Canonical Quantization 

We begin by highlighting the details of the canonical quan- 
tization of a scalar field pertinent to a perturbative treatment of 
non-equilibrium thermal field theory. 

The time-independent Schrodinger-picture field operator, 
denoted by a subscript S, may be written in the familiar plane- 
wave decomposition 

(1) 
where ^(p) - yp^ + M^ and flj(p; F,) and as(p; ti) are the usual 
single-particle creation and annihilation operators. It is essen- 
tial to emphasize that we define the Schrodinger-, Heisenberg- 
and Interaction (Dirac)-pictures to be coincident at the finite 
m/croscopic boundary time f,, i.e. 



<I)s(x;f,) = <l)H(f/,x;f,) = <l>i(f,-,x;f;) 



(2) 



It is at this picture-independent boundary time f, that initial con- 
ditions must be specified. The dependence upon the boundary 
time ti is separated from other arguments by a semi-colon. 

The time-dependent interaction-picture operator <E)i(x; f,) is 
obtained via the unitary transformation 



<Di(jc; ti) = e''^'* - ^>Os(x; f,)e-''^(-"= - ^> 



(3) 
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where H^^ is the free part of the Hamihonian in the Schrodinger 
picture. This yields 



J (Iny 2£(p) ^ 



(4) 



where ai(p, xq; If) - ai(p, 0; f,) g^'^^P)'" and its Hermitian conju- 
gate are the time-dependent interaction-picture annihilation and 
creation operators. These operators satisfy the canonical com- 
mutation relation 



[ai(p,.xo;f,), fl|(p',Xo;f,)] 



(5) 



with all other commutators vanishing. Note the presence of an 
overall phase e^'^'PX-to - 4) jjj pi fQj- j^^^ _£ ^^.^ 

In quantum statistical mechanics, we are interested in the 
Ensemble Expectation Values (EEVs) of operators at a fixed 
m/croscopic time of observation tj. Such EEVs are obtained by 
taking the trace with the density operator p(fy; f,), i.e. 



<•), = Z-\t)Txp{tf-ti)> 



(6) 



where Z,{t) - Tipitf', U) is the partition function, which is 
time-dependent in the presence of external sources. We have 
introduced the macroscopic time t - tf — U, which is the interval 
between the m/croscopic boundary and observation times. 

Consider the following observable, which is the EEV of a 
two-point product of field operators: 

0{x, y, tf- Id = Z-i (0 Trp(F/; Id^{If, x; f,)0(f^, y; f~) . (7) 

As shown in |[T1, it is not necessary to specify the picture in 
which the operators of the RHS of ([Til are to be interpreted, 
since all operators are evaluated at equal times. In addition, the 
observable O is invariant under simultaneous time translations 
of the boundary and observation times and depends only on the 
macroscopic time t: 0(x, y, Ij-; f,) = 0{x, y, If-I,', 0) = 0{x, y, t). 
Notice that O depends upon 7 independent coordinates: the spa- 
tial coordinates x and y and the macroscopic time t. 

The density operator p(f/; ?,) of a time-dependent and spa- 
tially inhomogeneous background is non-diagonal in the Fock 
space and contains an intractable incoherent sum of all possi- 
ble n to m multi-particle correlations, see El. We may account 
for our ignorance of the exact form of this density operator by 
defining the bilinear EEVs 

<fli(p,f/;f>|"(p'J/;f,)>, = (2;r)3 2£(p)<J«(p-p') 



+ 2£2(p)£2(p')/(p,p',0, (8a) 
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(alip'JfJdaiipJfM, = 2£2(p)£2(p')/(p,p',0, (8b) 

consistent with the canonical commutation relation (|5]l, where 
/(p, p',f) = /*(p',p, f)- The statistical distribution function 
/(p, p', t) is related to the particle number density n(q, X, t) via 
the Wigner transform 

«(q, X, t) = [^ e'Q-^ /(q + Q/2, q - Q/2, t) , (9) 



where we have introduced the relative and central momenta 
Q = p - p' and q = (p + p')/2, conjugate to the central and 
relative coordinates X = (x-i-y)/2 andR = x - y, respectively. 
Observe that spatial homogeneity is broken by the explicit de- 
pendence of /(p, p', t) on the two three-momenta p and p'. In 
the thermodynamic equilibrium limit, we have the correspon- 
dence/(p,p',0 ^ /eq(p,p') = (In)' 6^'\p - p')ME(p)), 
where fsix) - (e^'' - 1)"' is the Bose-Einstein distribution 
function andjS is the inverse thermodynamic temperature. 

3. Schwinger-Keldysh CTP Formalism 

We require a path-integral approach to generating EEVs for 
products of field operators. Such an approach is provided by 
the Schwinger-Keldysh CTP formahsm |47, 48 1. 

In order to obtain a generating functional of EEVs, we insert 
unitary evolution operators to the left and right of the density 
operator in the partition function 2(0 = Tr p(f/; If), yielding 

Z[p,J±,t] 

= Tr[te-'/",^'^^-«'^"WlpH(f/;f,)[Te''^"''''-*'^^^'*"^^T , (10) 

in the Heisenberg picture, where Q, is the temporally-bounded 
spacetime hypervolume [-f/2, f/2] x E^. We stress that ( 10 1 



differs fundamentally from existing interpretations of the CTP 
formalism |49, 50 1. Specifically, the Heisenberg -picture den- 
sity operator pH(f/; f,), which is explicitly time-dependent in the 
presence of the external sources J-t, is evaluated at the time of 
observation If and not the initial time f,-. In our approach, the 
role of the unitary evolution operators is to enable us to gen- 
erate EEVs for products of field operators as given in (JTl) by 
functional differentiation with respect to the external sources. 
The resulting EEVs are evaluated at the time of observation. 



We may interpret the evolution operators in ( 10 1 as defining 
a closed contour C = C+ U C- in the complex-time plane 
(t-plane, i € C), as shown in Figure [T] which is the union of 
two anti-parallel branches: C+, running from I, to If - ie/2; 
and C-, running from If - ie/2 back to f, - is. A small imag- 
inary part e = 0^ is added to separate the two, essentially co- 
incident, branches. We may introduce an explicit parametriza- 
tion of this contour z(u) [JJ, where u increases monotonically 
along C, which allows the definition of a path-ordering operator 
Tq. We emphasize that, in our formalism, this contour evolves 
in time, with each branch having length t. 

Following the notation of 1491 l50l . we denote fields con- 
fined to the positive and negative branches of the CTP contour 
by <l>±(x) = (I)(x" e C±,x). We then define the doublets 

(l)«(;c) = (o+(;t;), <l)_(x)), (11a) 

<1>„(X) = llab<i>\x) = ((t+(^), -<D-(X)), (lib) 

where the CTP indices a, b = 1,2 and riab - diag (1, -1) is 
an§0(l, 1) 'meti-ic' 



Inserting into ( 10 1 complete sets of eigenstates of the Heis- 
enberg field operator, we derive a path-integral representation 



Im t 



i(0) = U 
2(1) = ii - Je 



C_ 



Ret 

2(1/2) = if - ie/2 



macroscopic time t = 'R.ez{u) — ti 



initial conditions: 
macroscopic time i = 



observation: 
macroscopic time t = tf — ti 



Figure 1: The closed-time path, C - C+UC-. The relationship 
between microscopic and macroscopic times is indicated by a 
dashed black arrow. 



of the CTP generating functional HI, which depends on the 
path-ordered propagator 



iA"\x,yJf;ti) = {Tc[<i>"(x;td<I>''(yJd]), 

[ Af(x, y, if-, ti) A<(x, y, tf, h) 
' [A>(x, y. If, ti) Ad(x, y, tf ti) 



— I 



(12) 



For x^, y^ e C+, the path-ordering T^ is equivalent to the 
standard time-ordering T and we obtain the time-ordered Feyn- 
man propagator /AF(x,y, tf ti). On the other hand, for x", y° G 
C-, Tc is equivalent to anti-time-ordering T and we obtain the 
anti-time-ordered Dyson propagator /AD(x,y, f/;f,). For x° G 
C+ and y^ G C , x" is always 'earlier' than y^, yielding the 
absolutely-ordered negative-frequency Wightman propagator 
/A<(x,y, tf-, ti). Conversely, for y^ G C+ and x" G C-, we obtain 
the positive-frequency Wightman propagator /A>(x,y, tf ti). 

By means of a Legendre transform of the CTP generat- 
ing functional [IJ, we derive the respective Cornwall-Jackiw- 
Tomboulis effective action ifSTl . from which the CTP 
Schwinger-Dyson equation 



1 0,-1 



A„^'(x,3;,f/;f,) = A);^ '(x,}') + nat{x,y,tfti) 



(13) 



is obtained, where A^^'(x,y, Fy;f,) and A^^" (x,y) are the re- 
summed and free inverse CTP propagators, respectively, and 
nai,(x, y, tf ti) is the CTP self-energy, analogous in form to ( 12 1. 



4. Master Time Evolution Equations for 
Particle Number Densities 

In order to count both on-shell and off-shell contributions 
systematically, we 'measure' the number of charges, rather than 
quanta of energy. This avoids any need to identify 'single- 
particle' energies by means of a quasi-particle approximation. 
We begin by relating the Noether charge 



Q{xa;ti) = -/ I d^x(;rH(x;f,)<l>H(x;f;) - H. c.) 
to a charge density operator (3(q, X, Xq; ti) via 



(14) 



QiXa-ti) 



h'^m 



aq,X,Xo;f,). (15) 



By taking the equal-time EEV of (3(q, X,Xo; f,) and extracting 
the positive- and negative-frequency particle components, we 
arrive at the following definition of the particle number density 
in terms of off-shell propagators LJJ: 



"(-.-"•''^^-^/SB 



Q „-. 



QX 



X 0iqo)qoiA<iq+§,q-§,t;O) 



(16) 



using the translational invariance of the CTP contour. 

By partially inverting the CTP Schwinger-Dyson equation 



in ( 13 1, we derive the following master time evolution equation 
for the statistical distribution function /(q + j, q - j,t) |[T|: 

5,/(q+§,q-§,0 

- 2Jj^ ^ e-^^-q-Q0(qo)AJq+lq-lt;O) 



//g^^.-'^°'M^(^ 



+ ^,q-^,t;0) 



-If 



+ ^'(q-§,q+§,t;Q)) 

dqo dQo 

In 7.n 



?'?+§ 



'e°'%o)('^(? 
,f;0)). 



+ f,?-i,f;0) 



(17) 



where we have introduced 

.^(q+^,q- f,f;0) 



= -/; 



2 
d^fc 

(2^ 



4 iWriq +§,k,t;0) /A<(/fc, q-§,t;0) , 



(18a) 



'^(q+§,q-§,t;0) 



= lj^4 [in>{q+lk,t;Q)iAAk,q-§,t;Q) 
- iUJq +^,k,t;0)[ iAJk, q-^,t;0) 

- 2iAr{k,q-§,t;0))]. (18b) 

It is important to emphasize that ( [T7| provides a self-consistent 
time evolution equation for / valid to all orders in perturbation 
theory and to all orders in gradient expansion. The terms on the 



LHS of (17 1 may be associated with the total derivative in the 
phase space (X, p), which appears in the classical Boltzmann 



transport equation 1521. The expression ^ in ( 18a 1 is the force 
term, generated by the potential due to the dispersive part of the 



self-energy, and the "^ in ( 18b 1 are the collision terms. 



5. Non-Homogeneous Diagrammatics 

Let us consider a simple scalar theory, with one heavy real 
scalar field O and one light pair of complex scalar fields Ot' ' A')' 
described by the Lagrangian 

£ = i^^cD^O - \mW + d^^'d^x - mVx " 8<i>x''x -■■■ , 

(19) 
where the ellipsis contains omitted self-interactions. This model 
yields the following set of modified Feynman rules: 



Propagator 


Double-Momentum Representation 


Feynman (Dyson) 


+ 2n\2pQ\''l^6{p^ - M^)f(p, p', t)e'^P«-P'^'Vf2n\2p[j^l^6{p'^ - M^) 


+(-)ve-freq. Wightman 


'^>«)(/''/''' f/; ti) = 2n6(+(-)po)6{p^ - M^)(2nf5^\p - p') 

+ 2n\2po\^'^6{p^ - M^)f{p, p', t)e'^P''-P''^'>'f2n\2p[j^l^S{p'^ - M^) 


Retai-ded (Advanced) 


<^)^^'^'^N.o.(->e)^-p-M2^^^^''''^^-^'^ 


Pauli-Jordan 


i^\p,p') = 2ne{p^)6{p^ - M^){2nf6(^Hp - p') 


Hadamard 


ih\{p, p', tf- 1,) = 2nS{p^ - M^)(2nf6^^Hp - p') 

+ 27T\2po\^'^6{p^ - M^)2fip,p', t)e''^Po-PoMn\2p'g\^'^6(p'^ - M^) 


Principal-part 


iA''f,(p,p') = P ^_^^^(2;r)V'>(p-/) 



Table 1: The non-homogeneous free scalar propagators, wliere/(/:i,/?',f) = 9(po)9(p'Q)f(p,p',t) + 6{-po)9(-p'f^)f*(-p,-p\t). 



sum over all topologically distinct diagrams at a given 
order in perturbation theory. 

assign to each O-propagator line a factor of 

The set of non-homogeneous free propagators is listed in 
Table [U 

assign to each ;t'-propagator line a factor of 



= iA''-"''(p,p',tfJi) 



The double lines occurring in the CTP propagators re- 
flect the violation of three-momentum due to the non- 
homogeneous statistical distribution function /(p, p', f)- 

assign to each three-point vertex a factor of 



Pi 




-i8riatc(2nr6r(ZUPd 



where 77fl/,c... = I, a - b = ■■■ - I; rjahc- - -1, 
a = b = ■ ■ ■ - 2 and riahc- = otherwise. Due to the 
finite upper and lower bounds on all time integrals, the 
energy-momentum delta function is replaced by 

Sf\ ZL Pi) ^ S,i ZL PoM'\p - p') (20) 

in which energy conservation is systematically violated 
by the analytic weight function 

MZL Poj) ^ ^ sine [(2^1 po,,)f/2] . (21) 



This violation of energy conservation, shown diagram- 
matically by the dotted line terminated in a cross, results 
from the uncertainty principle, since the observation of 
the system is made over a finite time interval. 

associate with each external vertex a phase 

e'P^'f , 

where po is the energy flowing into the vertex. This phase 
results from the proper consideration of the Wick con- 
traction and field-particle duality relations. 

contract all internal CTP indices. 

integrate with the measure 

■dV 



J (2nf 



over the four-momentum associated with each contracted 
pair of CTP indices. 

• consider the combinatorial symmetry factors, where ap- 
propriate. 

These non-homogeneous Feynman rules encode the absolute 
spacetime dependence of the system starting from tree level. 

6. Absence of Pinch Singularities 

The perturbation series built from the non-homogeneous 
Feynman rules in Section l5] are free of the pinch singularities 
previously thought to spoil such perturbative treatments of non- 
equilibrium field theory, see e.g. 125] I53l455l . In our formula- 
tion, this absence of pinch singularities is ensured by two fac- 
tors: (i) the violation of energy conservation at early times and 




Figure 2: The one-loop €> self-energy !n^^^(^, q', tj; ti). 

(ii) the statistical distribution functions in free CTP propaga- 
tors are evaluated at the time of observation. The latter (ii) is 
in contrast to existing approaches in which free propagators do 
not evolve and depend only on the initial distributions. 

Consider the following one-loop insertion to the propagator: 

iA'-'^'"''(p,p',if;td = iA"'"''(p,p',tfJd 
+ iA^-^'ip, q. If, tdin^^](q, q'. If, IdiA^^'^iq', p' , If, ?,) . (22) 

Potential pinch singularities arise from terms like 



6(p^ - M^)6(po - Po)6(p'^ - M^) 



(23) 



However, at early times, energy is not conserved through the 
loop insertion. As a result, these terms are analytic, becoming 



dip' - M^)6>(po ■ 



■ p'o)S(p'^ 



■M^) 



(24) 



where S,{po - p[^) is given in (21 1. At late times, f — » co, 

Um 6,(po - p'o) = dipo - po) (25) 

and energy conservation is restored. However, in the same 
limit the system must have thermalized. In this case, the sta- 
tistical distribution functions appearing in free propagators will 
be the equilibrium distributions for which pinch singularities 
are known to cancel by virtue of the Kubo-Martin-Schwinger 
(KMS) relation 1561 . At intermediate times, pinch singulari- 
ties grow like a power law in t, which will always occur more 
slowly than the exponential approach to equilibrium. Thus, the 
perturbation series is free of pinch singularities for all times |[T]. 
Given the systematic diagrammatics of this approach, we 
may therefore truncate the master time evolution equations in 
( 17 1 in a perturbative loopwise sense. If the statistical distri- 
bution functions are tempered for all times, any ultra-violet di- 
vergences may be renormalized by the usual zero-temperature 
counter-terms, whilst infra-red divergences may be regularized 
by the partial resummation of thermal masses, see IJ\. 



7. Time-Dependent One-Loop Width 

To illustrate the distinctive features of our perturbative for- 
malism, let us consider two isolated but coincident subsystems 
S^(t, and ^, both separately in thermodynamic equilibrium and 
at the same temperature T = 10 GeV with the interactions 




1/ u/t 

X 
(b) 3 — » total annihilation 

X 




1/ u/t 

^2 \ / E2 

X X 
(c) 2 — > 1 Landau damping 

Figure 3: The four evanescent processes contributing to the 
one-loop time-dependent O width. 



switched off. The subsystem ,5^ij, contains only the field O 
with mass M - \ GeV and ,5^^, only the x fields of mass 
m - 0.01 GeV. At f = 0, we turn on the interactions and allow 
the system 5^ — 5^q, U ^ to re-thermalize. 

The one-loop non-local (b self-energy is shown in Figure l2j 
Neglecting back-reaction on the subsystem ^, the one-loop 
time-dependent O width is then given by the following integral: 



-d) 



r:,'(^,0 



647r3M 
X sinc[(g'o 



2 / 



d^k^^ 
E1E2 



Q'1,Q'2 = ±1 ' 

aiEi - a2E2) t] (1 + fsiaiEi) + fBiaiEz)) , 

(26) 



where £1 = ^^(k) = Vk^ + nf- and £2 = E^(q - k). The vi- 



olation of energy conservation, due to the sine function in ( [26) l, 

leads to otherwise-forbidden contributions from ai, a2 = -1 
(total annihilation) and a\ - -02 (Landau damping). In addi- 
tion, the kinematically-allowed phase space for 1 — > 2 decays is 
expanded. These evanescent processes are shown in Figure [3] 
where we have defined the evanescent action 



u = (qo- aiEi - aa^a) t , 



(27) 



quantifying the degree of energy non-conservation. Forf — > 00, 
we recover the known equilibrium result, since 



lim -sinc[(g'o-Q'i£'i-Q'2£'2)f] - 6{qo-aiEi-a2E2) . (28) 



In Figure |4j we plot the ratio 



r!-'(lql,0 



-d) 



rip'dqLO 



C(lql,r 



00) 



(29) 



of the time-dependent one-loop O width to its late-time equilib- 
rium value as a function of Mf for q^ - M^. In addition, we plot 
the separate contributions of the processes shown in Figure [3] 



\q\=10 GeV 



1.3 - 
1.2 - 

b; 1.1 - 
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0.9 - 

■ f 
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\q\ = 10 GeV 
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Figure 4: Left: the ratio P^' versus Mt for on-shell decays with |q| = 1 GeV (sohd black), 10 GeV (blue dotted) and 100 GeV (red 
dashed). Right: separate contributions to fjj,' for |q| = 10 GeV. Landau-damping contributions are equal up to numerical errors. 



8. Non-Markovian Oscillations 

In Figured we observe that the oscillations in the C) width 
have time-dependent frequencies. This non-Markovian beha- 
viour is inherent to truly out-of-equilibrium quantum systems, 
exhibiting so-called memory effects. Moreover, due to the 
Lorentz boost of ultra-violet modes relative to the rest frame of 
the heat bath, these memory effects persist for timescales much 
longer than the 1 /M which is to be expected for effects resulting 
from the uncertainty principle. 



In terms of the evanescent action m in ( 27 1 and in the high- 
temperature limit r » M, we may show quantitatively that the 
frequencies of these non-Markovian oscillations are given by 



„W 



(q, u, t) 



qo- 



(qlit) - Iqr +m]- ml)q,(t) + bag\q\A''^iql{t) - \q\\ m\, mj) 



2{ql(t) - IqP) 



(30a) 



0J2 \q,u,t) - 

iqlit) - IqP - ml + ml)qM - bae\q\A''\ql(t) - \q\\ m\, mj) 



2iql(t) - IqP) 



(30b) 



where fo, ag - ±1, A(x,y,z) = (x^ -y^ - z^)^ - 4-y^z^ andwe 

have introduced the evanescent energy 

u 
quit) = qo - - . (31) 

Notice that in the limit f — > 00, g'„(f) — > ^y and we ob- 
tain the usual time-independent kinematics. To the best of our 
knowledge, such a quantitative analysis of the non-Markovian 
evolution of memory effects has not been reported previously 
in the literature. 



9. Loopwise-Truncated Time Evolution Equations 



Truncating the master time evolution equation (17 1 to lead- 
ing order in a perturbative loopwise expansion, we obtain the 
following one-loop transport equation for the (^ statistical dis- 
tribution function: 



d,M\q\,t) = -- 



^ J (2n 



1 



1 



1 



■)3 2£a,(q) 2E^(k) 2E^{q - k) 



X — sinc[(Q'£'o(q) - aiE^ik) - a2E^(q - k))f/2J 

X (;r + 2Si[(Q'£'o(q) + a'i£'^(k) + aiE^iq - k))f/2]} 

X m-a) + Mqlt)] 

X [6(a,){\ + /^(|k|, 0) + e(-ai)f^(\k\, 0] 

X [0(a2)(l + f^(\q - k|, 0) + 0(-ff2)/^(lq - k|, t)] 

- [0(a)+Mq\,t)] 

X [0(o'i)/^(|k|, t) + e(-ai)(l + ./f (|k|, 0)] 

X [0(a2)f^(\q - k|, t) + 0(-a2)(l + /^(|q - k|, t))]} , 

(32) 



where a, U], 02 - ±1. The second and third lines of ( |32| i en- 
code the early-time violation of energy conservation. Replacing 
these lines by the Markovian approximation 



2;r6'(a)(5(£o(q) - d^Kk) - a2E2{q - k)) 



(33) 



we recover the semi-classical Boltzmann equation. However, 
given the equilibrium initial conditions of our model, this artifi- 
cial imposition of energy conservation along with the properties 
of the Bose-Einstein distribution ensure that the RHS of (|32li is 



zero for all times. Thus, the semi-classical Boltzmann equation 
cannot describe the re-thermalization of our simple model. This 
is true also for gradient expansions of Kadanoff-Baym equa- 
tions when truncated to zeroth order in time derivatives. Hence, 
it is only when energy-violating effects are systematically con- 
sidered, as in this new perturbative approach with all gradients 
included, that the dynamics of this re-thermalization is properly 
captured. 



It is clear that ( 32 1 describes only decay and inverse de- 
cay processes in the topologies shown in Figure l3] However, 
higher-multipUcity decays and scatterings can be systematically 
incorporated by consistently truncating the master time evolu- 



tion equation in ( 17 1 to a higher number of loops. 



10. Conclusions 

We have obtained master time evolution equations for par- 
ticle number densities that are valid to all orders in perturba- 
tion theory and to all orders in gradient expansion. The under- 
lying perturbation series are built from non-homogeneous free 
propagators and explicitly time-dependent vertices. Due to the 
systematic treatment of finite boundary and observation times, 
these diagrammatic series remain free of pinch singularities for 
all times. We are therefore able to truncate the time evolution 
equations in a perturbative loopwise sense, whilst keeping all 
orders in gradient expansion and capturing the dynamics on all 
timescales. This includes the prompt transient behaviour, which 
we have shown to be dominated by energy-violating processes 
that lead to non-Markovian evolution of memory effects. By 
virtue of our approach, we have been able to provide the first 
quantitative analysis of these memory effects. 

The foreseeable applications of this new formalism span 
high-energy physics, astro-particle physics, cosmology and con- 
densed matter physics. Dedicated studies of such applications 
will be the subject of future works. 
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